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Abstract 

This paper reviews results about free convection near a vertical flat plate em- 
bedded in some saturated porous medium. We focus on a third order autonomous 
differential equation that gives a special class of solutions called similarity solutions. 
Two cases are under consideration: in the first one we prescribe the temperature on 
the plate and in the second one we prescribe the heat flux on it. We will also see 
that the same equation appears in other industrial processes. 



1 Introduction 

Free convection boundary layer flows near a vertical flat plate embedded in some porous 
medium are studied for many years and a natural way to describe the convective flow is 
to look for similarity solutions. We consider two different sets of boundary conditions for 
the temperature on the plate: either we prescribe the temperature or we prescribe the 
heat flux. Both cases are leading to the same following third order non-linear autonomous 
differential equation 

f< + aff"-Pf* = (1.1) 

with the boundary conditions 

/(0) = -7, /' (00) = and /' (0) = 1, (1.2) 



or 

/(0) = -7, /' (00) = and /"(0) = -1. (1.3) 
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The first set of boundary conditions (jl.2|) with a = and (3 = m for m G R corresponds 
to prescribed heat on the plate as in [I], [5], [12], [E], [IHj, [21] and [21] • The second 
set of boundary conditions (jl.3j) with a = m + 2 and /3 = 2m + 1 for m G R is for the 
prescribed surface heat flux as done in [TOJ and [13]. In both cases the solutions depend 
on two parameters: m, the power-law exponent and 7, the mass transfer parameter. For 
7 = we have an impermeable wall, 7 < corresponds to a fluid suction, and 7 > to a 
fluid injection. 

Equation with suitable boundary conditions also arises in other industrial processes 
such as boundary layer flow adjacent to stretching walls (see [T], [2], [TH], [20], [22]) or 
excitation of liquid metals in a high-frequency magnetic field (see j2S])- 



2 The case of prescribed heat 
2.1 Derivation of the model 

We consider a vertical permeable flat plate embedded in a porous medium at the ambient 
temperature and a rectangular Cartesian co-ordinate system with the origin fixed at 
the leading edge of the vertical plate, the x-axis directed upward along the plate and the 
?/-axis normal to it. If we suppose that the porous medium is homogeneous and isotropic, 
that all the properties of the fluid and the porous medium are constants, that the fluid is 
incompressible and follows the Darcy-Boussinesq law and that the temperature along the 
plate is varying as x m the governing equations are 

du dv 
dx dy 

k f dp 

u = — 77- + pg 



/i \dx 

k dp 
v = 7T' 

dT dT x fd 2 T d 2 T 

U— + V— = X —r + 



dx dy \ dx 2 dy 2 
P — Ax>(l — P{T — Too)) 

where u and v are the Darcy velocities in the x and y directions, p, /1 and (3 are the 
density, viscosity and thermal expansion coefficient of the fluid, k is the permeability 
of the saturated porous medium and A its thermal diffusivity, p is the pressure, T the 
temperature and g the acceleration of the gravity. The subscript 00 is used for a value 
taken far from the plate. In our system of co-ordinates the boundary conditions along the 
plate are 

v(x,0) =ux m ^ 1 , T{x,0) =T w (x) =T OD + Ax m , m G R, 

with A > and ui G R (uj < corresponds to a fluid suction, u = is for an impermeable 
wall and u > corresponds to a fluid injection). The boundary conditions far from the 
plate are 

u(x, 00) = 0, T(x, 00) = Too. 
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If we introduce the stream function ip such that 

dip dip 
dy ' dx 

and assuming that convection takes place in a thin layer around the heating plate, we 
obtain the boundary layer approximation 



= Poo(3gk dT 
dy 2 /j, dy ' 

d 2 T 1 fdTdip dTdip 



dy 2 X \dx dy dy dx 
with 

dip m-i dip 

— — [x, 0) = — cox 2 and — — [x, oo) = 0. 

dx dy 

Let us introduce the new dimensionless similarity variables 

t=(Ra x )^, iP(x,y) = \(Ra x ) l *f(t), 6{t) /( ' r '' /) ~- / - 



(2.1) 



(2.2) 



with i?a x = (pooPgk(T w (x) — T oc )x)/( / uA) the local Rayleigh number. In terms of these 
variables equations (|2.1j) and ()2.2j) become 

/" - 6' = 0, (2.3) 

and 

(9" + - mf'O = 0, 

with the boundary conditions 

/(0) = - T , 0(0) = 1, 

and 

/'(oo) = 0, 0(oc) = O, (2.4) 
where the prime denotes differentiation with respect to t and 



2a; / IT 
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m + 1 y poofigkAX 

Integrating ()2.3|) and taking into account the boundary conditions ()2.4|) leads to 

f = e 

and the problem (|l.lj) - (|1.2|) with a = and j3 = m for m G K follows. 
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2.2 Usefull tools 



2.2.1 The initial value problem 

Let -P m , 7i a be the following initial value problem 

/"' + ^ff" - mf' 2 = 

/(0) = ~ 7 ' (2 5) 

/' (o) = i, [b) 

f{0) = a. 

This first approach used is a shooting method that consists in finding values of f"(0) = a 
for which / exists on [0, oo) and such that /'(oo) = 0. This direct method allows us to 
consider vanishing solutions but does fail in some cases (see jHj). 



2.2.2 The blowing-up co-ordinates 

Let us notice that if / is a solution of then for all k > the function t i — > Kf(Kt) 
is a solution too. Then, considering an interval I on which a solution / of (jl.lj) does not 
vanish, for r G I we can introduce the following blowing-up co-ordinates 

We/, 8 = J*f{£)d£, »W = Jffi ^d M = J$- ( 2 -6) 
Then, we easily get 

J u = P(u, v) := v - 2m 2 , , . 

1 u = Q m (u,v) := -^ + mti 2 -3w, 1 J 

where the dot is for differentiating with respect to the variable s. To come back to the 
original problem it is sufficient to consider the initial value problem P m ,-y,a with 7 ^ and 
look at the trajectories of the corresponding plane dynamical system ([2.7)1 . For details, 
see [TTj. 



2.3 Main results 

The problem (|l.l|) - (jl.2J) appears in engineering and physical litterature, in very different 
context, in the middle of the previous century. 

Rigorous mathematical results arise around the sixties. In |2E] (Appendix 2) it is men- 
tionned that for 7 = a simple explicit solution can be obtained in both cases m = 1 and 
m — — I (see also [21], jl] and [3]). On the other hand, the author notes that Mr J. 
Watson has given a simple proof that (jl.l|) - (jl.2|) has no solution for 7 = and m < —1. 
An explicit solution is also given for m = 1 and any 7, first in [2H], and later in j2Z| and 
[Hj. In these latter papers one can also find the explicit solution for m — — | and any 7. 
Nonexistence for 7 = and m = — | was noted in jl]. In £2], it is shown that for 7 = 
and m < — | there are no solutions satisfying /'/ 2 — > at infinity. 

Recently, further mathematical results concerning existence, nonexistence, uniqueness, 
nonuniqueness and asymptotic behaviour, are obtained in jS] for 7 = 0, and in |T8] . 
[Hj, \H\ [IH] and [0] for the general case. 
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Numerical investigations can be found in [T], 0, [12], [E], [2J], [22] and [2Zj- 

In view of all these papers, the following conclusions can be drawn. 



• For m < — 1, there exists 7* > such that problem (jl.lj) - (jl.2|) has infinitely many 
solutions if 7 > 7*, one and only one solution if 7 = 7*, and no solution if 7 < 7*. 
For 7 = 7* we have that f(t) — > A < as t — > 00, and for every 7 > 7* there are 
two solutions / such that f(t) — > A < as £ — > 00 and all the other solutions verify 
/(i) — > as t — > 00 Moreover, if / is a solution to (|l.lj) - (|1.2|) . then / is negative, 
strictly increasing and either concave or convex-concave. (See Fig. 1 for the two 
solutions such that f(t) — > A < as t — > 00 and three other solutions in the case 
m = —2 and 7 = 5.) 

-2.5 
-3 
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-4 
-4.5 

-5 

5 10 15 20 

Fig. 1 

• For m = — 1 and for every 7 e R, the problem (jl.lj) - ()1.2|) has no solution. 

• For —1 < m < — I and for every 7 > 0, the problem (jl.lj) - (jl.2|) has no solution. 

• For — 1 < m < — |, there exists 7* < such that problem (|l.lj) - ([1.2)) has no solution 
for 7# < 7 < 0, one and only one solution which is bounded for 7 = 7*, and 
two bounded solutions and infinitely many unbounded solutions for 7 < 7*. These 
solutions are strictly increasing and either concave or convex-concave. (See Fig. 2 to 
see the two bounded solutions and four unbouded solutions in the case m = —0.75 
and 7 = —10.) 

• For — I < m < — I and for every 7 < 0, the problem (jl.lj) - ()1.2j) has one bounded 
solution and infinitely many unbounded solutions. All these solutions are strictly 
increasing and either concave or convex-concave. 

• For — i < m < and for every 7 G R, the problem ()l.lj) - ()1.2JI has an infinite number 
of solutions. Moreover, if 7 < one and only one solution is bounded, and if 7 > 
at least one is bounded, many infinitely are unbounded. All solutions are strictly 
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increasing and either concave or convex-concave. If 7 > 0, the solutions becomes 
positive for large t. 
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2 4 6 8 10 

Fig. 2 

• For m G [0, 1] and for every 7 e R, the problem (jl.lj) - (J1.2j) has one and only one 
solution, moreover this solution is concave and bounded. (See Fig. 3 for the unique 
solution in the case m = 0.5 and 7 = 0.) 

1.4 i 1 1 1 1 1 

1.2 - 
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Fig. 3 

• For m > 1 and for every 7 G 1, the problem (jl.l|l - (|1.2jl has one and only one 
concave solution and an infinite number of concave-convex solutions. All these 
solutions are bounded. Moreover, there is an unique concave-convex solution that 
verifies f(t) — > A > as t — > 00 and all the other concave-convex solutions are such 
that f(t) —>■ as t —>■ 00. (See Fig. 4 for the unique concave solution and three 
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concave-convex solutions in the case m = 1.1 and 7 = 0.) 




2 4 6 8 10 



Fig. 4 

Remark 2.1 The case m = leads to the well know Blasius equation (see jTjj that 
also is a special case of the Falkner-Skan equation (see |17j ). 

Remark 2.2 In j2B| the authors gives results about a slightly different problem for m = 
— 1 that involves pseudo-similarity . 

We see from these results, that the unsolved questions concern the case 7 > 0. More 
precisely, it should be interesting to try to answer to the following points 

• For — ~ < m < — |, what happens for 7 > ? 

• For — I < m < and 7 > 0, is there one or more bounded solutions ? 
Another purpose is to compute the critical values 7* appearing in the results above. 

3 The case of prescribed heat flux 

We now suppose that the plate is subjected to a variable heat flux varying as x m and a 
mass transfer rate varying as x~ ~ following pHJ to obtain the problem (jl.lj) - (jl.3j) with 
a = m + 2 and (3 = 2m + 1. The mathematical study is made in ^0] and leads to the 
following results 

• For m < — 2 there exists 7* > ^/ ( m ^ 2 ) 2 sucn ^ na ^ ^ ne P r °bl em (jl.l|) - (|1.3J) has no 
solution for 7 < 7*, one and only one solution for 7 = 7* and infinitely many 
solutions for 7 > 7*. For 7 = 7* we have that f(t) — > A < as t — >• 00, and for every 
7 > 7* there are two solutions / such that f(t) — > A < as t — > 00 and all the other 
solutions verify f(t) — >• as t — >• 00 Moreover, if / is a solution of ljl.l)l - (|1.3|L then 
/ is negative, strictly concave and increasing. 
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• For m = —2 and for every 7 G R, the problem (|l.l)) - f|1.3jl has no solution. 

• For —2 < m < — 1, there exists 7* < such that the problem p.lj) - (jl.3j) has no 
solution for 7 > 7*, one and only one solution which is bounded for 7 = 7* and two 
bounded solutions and infinitely many unbounded solutions for 7 < 7*. Moreover, 
if / is a solution of (|l.ip - ([1.3|) . then / is positive, strictly concave, increasing and 

/ (0) > — ( m +2) 7 • 

• For m = — 1 the problem (|1.1)) - ([1.3|) only admits solutions for 7 < 0. In this case 
there is an unique bounded solution with /'(0) = — - and an infinite number of 
unbounded solutions with /'(0) > —-. Moreover all the solutions are positive, 
strictly concave and increasing. 

• For —1 < m < — | the problem (jl.l|) - (jl.3j) admits at least one bounded solution for 
7 G R and many infinitely unbounded solutions for 7 < 0. All these solutions are 
increasing and strictly concave and uniqueness of the bounded solution holds for 
7 < 0. 

• For m > — ~ all the solutions are bounded. 

• For — 2 < m < 1 and for every 7 G R the problem (jl.lj) - (jl.3|1 has one and only one 
solution. This solution is strictly concave and increasing. (Let us notice that for 
m = — \ we have the Blasius equation.) 

• For m > 1 and 7 G R the problem has one and only one concave solution 
and infinitely many concave-convex solutions. Moreover, there is an unique concave- 
convex solution that verifies f(t) — > A > as t — ► 00 and all the other concave- 
convex solutions are such that f(t) — > as t — > 00. 

In this case it remains only two open questions 

• For —1 < m < — ~ and 7 > 0, is the bounded solution unique ? 

• For —1 < m < — \ and 7 > 0, is there unbounded solution ? 

4 Asymptotic behaviour of the unbounded solutions 

For the equation (jl.lj) we have the following asymptotic equivalent found in |Hj and ^Hj 
that holds for unbounded solutions 

Theorem 4.1 Let f be an unbounded solution of or There exists 

a constant c > such that 

a 

\f(t)\ ~ cfr^P as t — > 00. 
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